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STATISTICS ON GRAPHS, EXPONENTIAL FORMULA AND 

COMBINATORIAL PHYSICS 

Laurent Poinsot, Gerard H. E. Duchamp, Silvia Goodenough and Karol A. Penson ™ 



Abstract. The concern of this paper is a famous combinatorial 
formula known under the name "exponential formula" . It occurs 
quite naturally in many contexts (physics, mathematics, computer 
science). Roughly speaking, it expresses that the exponential gen- 
erating function of a whole structure is equal to the exponential 
of those of connected substructures. Keeping this descriptive state- 
ment as a guideline, we develop a general framework to handle many 
different situations in which the exponential formula can be applied. 
Keywords. Combinatorial physics, Exponential generating func- 
tion, Partial semigroup, Experimental mathematics. 



1 Introduction 

Applying the exponential paradigm one can feel some- 
times incomfortable wondering whether "one has the 
right" to do so (as for coloured structures, for exam- 
ple). The following paper is aimed at giving a rather 
large framework where this formula holds. 

Exponential formula can be traced back to works by 
Touchard and Ridell & Uhlenbeck [2U1 [17]. For an other 
exposition, see for example [4j [7J [Qj [19] . 

We are interested to compute various examples of EGF 
for combinatorial objects having (a finite set of) nodes 
(i.e. their set-theoretical support) so we use as central 
concept the mapping a which associates to every struc- 
ture, its set of (labels of its) nodes. 

We need to draw what could be called "square-free de- 
composable objects" (SFD). This version is suited to our 
needs for the "exponential formula" and it is sufficiently 
general to contain, as a particular case, the case of mul- 
tivariate series. 
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2 Partial semigroups 

Let us call partial semigroup a semigroup with a partially 
defined associative law (see for instance [6] for usual semi- 
groups and [TJ HH [Tg] for more details on structures with 
a partially defined binary operation). More precisely, a 
partial semigroup is a pair (S, *) where S is a set and * 
is a (partially defined) function S X S — )• S such that the 
two (again partially defined) functions S X S x S — >• S 

(x, y,z) h-t (x * y) * z and (x, y , z) M- x * (y * z) (1) 

coincide (same domain and values). Using this require- 
ment one can see that the values of the (partially defined) 
functions S n — > S 

(xi, ■ ■ ■ ,x n ) H> E T (xi, ■ ■ ■ ,x n ) (2) 

obtained by evaluating the expression formed by labelling 
by Xi (from left to right) the ith leaf of a binary tree T 
with n nodes and by * its internal nodes, is independant 
of T . We will denote x\ * • • • * x n their common value. 
In this paper we restrict our attention to commutative 
semigroups. By this we mean that the value X\ * ■ • ■ * 
x n does not depend on the relative order of the a;,-. A 
nonempty partial semigroup (S, *) has a (two-sided and 
total) unit e € S if, and only if, for every u> £ S, uj * e — 
uj = e*oj. Using associativity of *, it can be easily checked 
that if S has a unit, then it is unique. 

Example 2.1. Let F be a set of sets (resp. which con- 
tains as an element) and which is closed under the dis- 
joint sum U, i.e., if A, B G F such that A n B = 0, then 
A U B(= A U B) e F. Then (F, U) is a partial semigroup 
(resp. partial semigroup with unit). 

3 Square-free decomposable par- 
tial semigroups 

Let 2^ N+ ) be the set of all finite subsets of the positive 
integers N + and (S, ©) be a partial semigroup with unit 
(here denoted e) equipped with a mapping a : S — > \ 
called the (set-theoretic) support mapping. Let D be the 
domain of ©. The triple (S, ©, <j) is called square-free 
decomposable (SFD) if, and only if, it fulfills the two fol- 
lowing conditions. 



2 



L. Poinsot, G. H. E. Duchamp, S. Goodenough and K.A. Penson 



• Direct sum (DS): 

1. <t(uj) = iff oj = e; 

2. D = {(w 1 ,w 2 )eS' 2 :CT(wi)n ( 7(cj2) = 0}; 

3. For all wi,W2 G 5, if (wj.,^) € -D then ct(wi © 
w 2 ) = cr(wi) U cr(w2)- 

• Levi's property (LP): For every oji,u)2,oj ,ur € 5 
such that (wi,w 2 ), (w 1 ,^ 2 ) G £> and U\ (Boj 2 = oj 1 ® 
oj 2 , there are G S for t = 1, 2, j = 1,2 such that 
(w*,a; 2 ), Wj) € D, oj 1 = Wj©w 2 and = a^ffiu^ 
for t = l,2 and j = 1, 2. 

Remark 3.1. The second and third conditions of (DS) 
imply that ct(oji®oj2) = a(oJi)Ua(u>2) whenever (wi,^) S 
D (which means that a(oJi)na(oj 2 ) = 0), where U denotes 
the disjoint sum. 

Example 3.1. As example of this setting we have: 

1. The positive square- free integers, cx(n) being the set 
of primes which divide n, the atoms being the prime 
numbers. 

2. All the positive integeres (S — N + ), under the usual 
integer multiplication, er(n) being the set of primes 
which divide n. 

3. Graphs, hypergraphs, (finitely) coloured, weighted 
graphs, with nodes in N + , a(G) being the set of nodes 
and © the juxtaposition (direct sum) when the set 
of nodes are mutually disjoint. 

4. The set of endofunctions / : F — > F where F is a 
finite subset of N + . 

5. The (multivariate) polynomials in N[X], X — {xi : 
i€/}, with I C N + , being a nonempty set of (com- 
muting or not) variables, with a(P) — Alph(P) the 
set of indices of variables that occur in a polynomial 
P, and © = +. 

6. For a given finite or denumerable field, the set of irre- 
ducible monic polynomials is denumerable. Arrange 
them in a sequence (Pn)neN+i then the square-free 
monic (for a given order on the variables) polynomi- 
als is SFD, a(P) := {n 6 N+ : P n divides P} and © 
being the multiplication. 

7. Rational complex algebraic curves; cr(V) being the 
set of monic irreducible bivariate polynomials van- 
ishing on V. 

In what follows we write ffi" =1 Wi instead of ui\ ® • • • @bJ n 
(if n = 0, then ©■L 1 u>i = e) and we suppose that (S, ffi, a) 
is SFD for the two following lemmas. 

Lemma 3.1. Let u>i, . . . ,ui n £ S such that ffi™ =1 o;.; is 
defined. Then for every i,j £ {1, . . . , n} such that i ^ j , 
it holds that o~(u>i) PI cr(wj) — 0- In particular, if none u>k 



is equal to e, then Wi ^ ojj for every i,j G {1, . . . , n} such 

n 

that i 5^ j. Moreover er(ffi™ =1 Wi) = | | <r(uji). 

»=l 

Lemma 3.2. Let (wi)™ =1 be a finite family of elements 
of S with pairwise disjoint supports. Suppose that for 
i = 1, • • • ,n, uJi = ®^.L 1 u>^, where (cj| : )^1 1 is a finite 
family of elements of S . Then ffi" =1 Wi = ©™ = i {^Li^i )• 

These lemmas are useful to define the sum of two or 
more elements of S using respective sum decompositions. 

Now, an atom in a partial semigroup with unit S is any 
object oj 7^ e which cannot be split, formally 

oj = oji © oj 2 e € {oji, oj 2 } ■ (3) 

The set of all atoms is denoted by atoms(5). Whenever 
the square-free decomposable semigroup S is not trivial, 
i.e., reduced to {e}, atoms(5) is not empty. 

Example 3.2. The atoms obtained from examples 13.11 

1. The atoms of 13.1121 are the primes. 

2. The atoms of 13.1131 are connected graphs. 

3. The atoms of 13.1141 are the endofunctions for which 
the domain is a singleton. 

4. The atoms of 13.1151 are the monomials. 

The prescriptions (DS,LP) imply that decomposition 
of objects into atoms always exists and is unique. 

Proposition 3.1. Let (S,®,a) be SFD. For each oj e 
S there is one and only one finite set of atoms A — 
{oji, ■ ■ ■ ,oj n } such that oj = (B™ =1 oji. One has A = 
iffu = e. 

4 Exponential formula 

In this section we consider (S, ®,cr) as a square- free 
decomposable partial semigroup with unit. 

In the set S, objects are conceived to be "measured" by 
different parameters (data in statistical language). So, to 
get a general purpose tool, we suppose that the statistics 
takes its values in a (unitary) ring R of characteristic zero 
that is to say which contains Q (as, to write exponential 
generating series it is convenient to have at hand the frac- 
tions -if). Let then c : S — > R be the given statistics. For 
F a finite set and each X C S, we define 

X F := {oj e X : u{u) = F} . (4) 

In order to write generating series, we need 

1. that the sums c(Xp) :— c ( a; ) exist for every 

finite set F of N + and every ICS; 
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2. that F — > c{Xp) would depend only of the cardinal- 
ity of the finite set F of N + , for each fixed ICS; 

3. that c(u)\ © u>2) = c(uj\) .c(oj2) ■ 
We formalize it in 

(LF) Local finiteness. — For each finite set F of N + , the 
subset Sp of S is a finite set. 
(Eq) Equivariance. — 

card(Fi) = card(F 2 ) ==> c(atoms(iS)j? 1 ) = c(atoms(5)jr 2 ) . 

(5) 

(Mu) Multiplicativity . — 

c(wi © W2) = c(wi).c(w 2 ) • (6) 

Remark 4.1. a) In fact, (LF) is a property of the set 
S, while (Eq) is a property of the statistics. In practice, 
we choose S which is locally finite and choose equivariant 
statistics for instance 

£^^j) ^(number of cyclcs)^(numbcr of fixed points) 

for some variables x, y. 

b) More generally, it is typical to take integer-valued 
partial (additive) statistics c±, ■ ■ ■ Ci, ■ ■ ■ , c r (for every 
u G S, Ci(w) G N) and set c(ui) = x? {u) x% {u3) ■ ■ ■ x c r A " ] . 

c) The set of example 13.1121 is not locally finite, but 
other examples satisfy (LF): for instance 13. 1131 if one asks 
that the number of arrows and weight is finite, 13.1111 

A multiplicative statistics is called proper if c(e) 7^ 0. 
It is called improper if c(e) = 0. In this case, for every 
ui G S, c(ui) = as c(lu) = c(uj © e) = c(w)c(e) = 0. 

If R is a integral domain and if c is proper, 
then c(e) = 1 because c(e) = c(e © e) = c(e) 2 , 
therefore 1 = c(e). Note that for each X C S, 

c(X$) = c ( w ) = { t llx- For every 

finite subset X of S, we also define c(X) := c(uS), 

LJ£X 

then we have in particular c(0) = (which is not the 
same as c(S$) — c({e}) if c is proper). The requirement 
(LF) implies that for every X C S and every finite set 
F of N + , c(Xp) is defined as a sum of a finite num- 
ber of terms because Xp C SV, and therefore Xp is finite. 

Now, we are in position to state the exponential formula 
as it will be used throughout the paper. Let us recall the 
usual exponential formula for formal power series in i2[[z]] 
(see [T3l [19] for more details on formal power series) . Let 
f(z) — >J fnZ n . Then we have 

n>l 

..n 

e f = J2 a A (7) 

n>0 



where 

a n = J, IT /card(p) (8) 

with LI„ being the set of all partitions of [l..n] (in partic- 

ular for n = 0, ao = 1) and e z = — - G 

n>0 

In what follows [l..n] denotes the interval {j G N + : 1 < 
j < n}, reduced to when n = 0. Let (5, ffi,cr) be a 
locally finite SFD and c be a multiplicative equivariant 
statistics. For every subset X of S one sets the following 
exponential generating series 

oo n 

EGF(I;z) = ^c(I M )^. (9) 

71=0 

thm 4.1 (exponential formula). Let S be a locally finite 
SFD and c be a multiplicative equivariant statistics. We 
have 

EGF(S*; Z) = C(e) - 1 + e EGF(atoms(S); Z ) ( 1Q ) 

In particular if c(e) = 1 (for instance if c is proper and 
R is an integral domain), 

EGF(S; Z) = e EGF(atoms(S);z) _ ^ 

Proof — Let n — 0. Then the unique element of S$ 
is e. Therefore c(5o) = c(e). Now suppose that n > 
and let cj G >Srx.. n ]. According to proposition 13.11 there 
is a unique finite set {ai, . . . , ctk\ Q atoms(S') such that 
w = ffi^ =1 a,. By lemma I3TT1 {cr(aj) : 1 < i < k} is a par- 
tition of [l..n] into k blocks. Therefore ui G atoms(S)p © 
• ■ • © atoms(S , ) Pfc where Pj = cr(aj) for i = 1, . . . , k. We 
can remark that a± © • • • © a.^ is well-defined for each 
(ai, . . . , Qife) G atoms(iS')p x • ■ ■ x atoms(S f ) Pfc since the 
supports are disjoint. Now, one has, thanks to the parti- 
tions of [l..n] 

S[i..n] = U 0atoms(5) p (12) 
c(5 [1 .. n] )= J] J] c(atoms(5) p ) (13) 

as, for disjoint (finite) sets F and G of N + , it is easy to 
check that c(X F © X G ) = c(X F )c(X G ) for every ICS 
and because the disjoint union as only a finite number of 
factors. Therefore due to equivariance of c on sets of the 
form atoms(S') F , one has 

C(S[1..„])= II C ( at0mS ( 5 )[l"Card(p)]) ' ( 14 ) 

•Ten,, pe-rr 

But c(atoms(S')[ 1 carc |(p)]) is the card(p)th coefficient of 
the series EGF(atoms(S'); z). Therefore due to the 
usual exponential formula, EGF(S';z) = c(e) — 1 + 

e EGF(atoms(S); Z )^ Nqw jf ^ = ^ then we obtain 
EGF(S'; Z) = e EGF(atoms(5);z)^ 
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5 Two examples 

The examples provided here pertain to the class of la- 
belled graphs where the "classic" exponential formula ap- 
plies, namely Burnside's Classed Burn a ^, defined, for 
< a < b two integers, as the class of graphs of numeric 
endofunctions / such that 

r = t (is) 

where /" denotes the nth power with respect to func- 
tional composition. Despite of its simplicity, there are 
still (enumerative combinatorial) open problems for this 
class and only Bii + i gives rise to an elegant formula 
[51 IT5] (see also [TT], for the idempotent case: £ = 1 and 
compare to exact but non-easily tractable formulas in [4] 
for the general case in the symmetric semigroup, and in 
[12] for their generalization to the wreath product of the 
symmetric semigroup and a finite group). 

The second example: the class of finite parti- 
tions which can be (and should here) identified as 
graphs of equivalence relations on finite subsets F C 
N + . Call this class "Stirling class" as the num- 
ber of such graphs with support [l..n] and k con- 
nected components is exactly the Stirling number of 
the second kind 52 (n, k) and, using the statistics 

^.(number of points) y(number of connected components) Qj^g ob- 
tains 

£ S 2 (n,k)^-y k = e y^ . (16) 

TV. 

n,k>0 

Examples of this kind bring us to the conclusion that 
bivariate statistics like Burn a ^{n, k), S2 (n, k) or 5i(n, k) 
(Stirling numbers of the second and first kind) are better 
understood through the notion of one-parameter group, 
conversely such groups naturally arinsing in Combinato- 
rial Physics lead to such statistics and new ones some of 
which can be interpreted combinatorially. 

6 Generalized Stirling numbers in 
Combinatorial Physics 

In Quantum Mechanics, many tools boil down to the con- 
sideration of creation and annihilation operators which 
will be here denoted respectively a) and a. These two 
symbols do not commute and are subject to the unique 
relation 

[0,0+] = 1. (17) 

The complex algebra generated by these two symbols and 
this unique relation, the Heisenberg-Weyl algebra, will be 
here denoted by HWc- The consideration of evolution 

1 The name is related to the notion of free Burnside semigroups, 
namely the quotient of the free semigroup A + , where A is a finite 
alphabet, by the the smallest congruence that contains the relators 
uj n + m = u) n , ui £ A+. For more details see 1151 . 



(one-parameter) groups e Af! where f2 = Oi(u>)u) is 

uieHWc 

an element of HWc, with all - but a finite number of 
them - the complex numbers a{oj) equal to 0, and to a 
word on the alphabet {a,a^} leads to the necessity of 
solving the Normal Ordering Problem, i.e., the reduction 
of the powers of SI to the form 

n n = 5^ J 8 i j(a t )V . (18) 

In the sequel, Normal{Vl n ) denotes such a sum. This 
problem can be performed with three indices in general 
and two in the case of homogeneous operators that is 
operators for which the "excess" e = i — j is constant 
along the monomials (a 1 ) 1 ^ of the support (for which 
foj ^ 0). Thus, for 

fi = /M at ) v (19) 

i— j=e 

one has, for all «eM, 

00 

Normal(n n ) = {o)) ne ^ 5 n (n, k){o)) k a k (20) 

k=0 

when e > 0, and 

Normal(n n ) = t ^ S n (n, k){a)) k a k J a n ^ (21) 
\fe=0 / 

otherwise. It turns out that, when there is only one an- 
nihilation, one gets a formula of the type (x, y are formal 
commutative variables) 

V Sn(n,k)^-y k = g(x)e^^ s ^ n ^ ( 22 ) 
* — ' nl 

n,k>0 

which is a generalization of formula (I16|) . A complete 
study of such a procedure and the details to perform the 
solution of the normal ordering problem may be found in 

0- 

7 Conclusion 

In this paper, we have broadened E| H the domain of appli- 
cation of the exponential formula, a tool originated from 
statistical physics. This broadening reveals us, together 
with the essence of "why this formula works" , a possibility 
of extension to denominators other than the factorial and, 
on the other hand, provides a link with one-parameter 
groups whose infinitesimal generators are (formal) vector 
fields on the line. The general combinatorial theory of 
the correspondence (vector fields «-» bivariate statistics) 
is still to be done despite the fact that we have already a 
wealth of results in this direction. 

2 A part of our setting can be reformulated in the categorical 
context [2j [3] 

3 Another direction is the q-exponential formula |10II16| . 
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